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Abstract A measurement of event shape variables is pre- 
sented for large momentum transfer proton-proton collisions 
using the ATLAS detector at the Large Hadron Collider. Six 
event shape variables calculated using hadronic jets are stud- 
ied in inclusive multi-jet events in 35 pb~ 1 of integrated lu- 
minosity at a center-of-mass energy of y/s = 7 TeV. These 
measurements are compared to predictions by three Monte 
Carlo event generators containing leading-logarithmic par- 
ton showers matched to leading order matrix elements for 
2—^2 and 2 — >• n (n = 2, ...6) scattering. Measurements of 
the third-jet resolution parameter, aplanarity, thrust, spheric- 
ity, and transverse sphericity are generally well described. 
The mean value of each event shape variable is evaluated as 
a function of the average momentum of the two leading jets 
Pt,\ and pr,2, with a mean pj approaching 1 TeV. 
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1 Introduction 

Event shapes represent a generic class of observables that 
describe the patterns, correlations, and origins of the en- 
ergy flow in an interaction. In terms of hadronic jet produc- 
tion, event shapes are an indirect probe of multi-jet topolo- 
gies. These observables have had a long and fruitful his- 
tory, having been used to measure the strong coupling con- 
stant as and to test asymptotic freedom [1^1], to constrain 
color factors for quark and gluon couplings [5], to assess 
the accuracy of leading order (LO) and next-to-leading or- 
der (NLO) Monte Carlo (MC) generators [6, 7], to deter- 
mine the contribution of non-perturbative quantum chromo- 
dynamics (QCD) power corrections [8], and to search for 
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physics beyond the Standard Model [9]. Furthermore, re- 
cent efforts to provide advanced, high-precision theoretical 
calculations of a range of event shapes for the Large Hadron 
Collider [10, 11] provide renewed impetus for making such 
measurements. 

This analysis considers six event shapes calculated using 
hadronic jets. These observables are crucially tied to both 
the multi-jet nature of the final state produced in high en- 
ergy collisions and have a strong history in the literature: the 
third-jet resolution parameter [4, 12-14], j23i the sphericity 
and transverse sphericity [15, 16], S and S±; the aplanarity, 
A; and the event thrust and its minor component [17], T± 
and T m x- Events with high transverse momentum central 
leading-jet pairs are used for the measurements. Each event 
shape variable is defined such that it vanishes in the limit 
of a pure 2 — > 2 process and increases to a maximum for 
uniformly distributed energy within a multi-jet event. Hard 
gluon emission is thereby signified by large non-zero values 
of each observable. Furthermore, some of the event shape 
variables are evaluated as ratios of final state observables, 
which reduces their sensitivity to jet energy scale (JES) cal- 
ibration uncertainties as well as other experimental and the- 
oretical uncertainties. These measurements permit detailed 
tests of the phenomenological models of QCD in leading 
order MC programs and to indirectly test the running of as 
through measurements performed as a function of the aver- 
age leading jet momentum. In addition, these results may be 
used to provide input to tune MC generators in the future. 

All event shapes measured in this analysis are defined 
using jets to represent the final state four-momenta, as dis- 
cussed in Section 2. The ATLAS detector is described in 
Section 3, with a particular emphasis on the components 
relevant for the measurement of event shape variables. Sec- 
tion 4 presents the event selection and description of simu- 
lated events which are compared to the data. Jet definitions, 
calibrations, and selection criteria are also described in Sec- 
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tion 4. Finally, the results of these measurements are pre- 
sented in Section 5. 



2 Event Shape Definitions 

Six event shapes are measured using high transverse mo- 
mentum (/;>t) jets. The first observable, V23, is a measure of 
the third-jet pj relative to the summed transverse momenta 
of the two leading jets in a multi-jet event and is defined as: 



V23 



Pt,3 
n T2 



(1) 



where Hj2 = (pr,i + Pt,i) is the scalar sum of jet momenta 
and the subscript i = 1,2,3 refers to the leading, subleading, 
or third leading jet in the event. The range of allowed values 
fory23 is < V23 < 1/4 and it is often expressed as lny23 [1 1, 
16]. This definition is different from the original [12] defini- 
tion which uses the JADE jet algorithm [18]. Eq. 1 is defined 
with an explicit third-jet as opposed to a continuously vari- 
able threshold in the jet algorithm. The sphericity, S, trans- 
verse sphericity, S±, and aplanarity, A, embody more global 
information about the full momentum tensor of the event, 
Mxyz, via its eigenvalues Ai, A2 and A3: 
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where the sum runs over all jets used in the measurement. 
The individual eigenvalues are normalized and ordered such 
that X\ > X2 > A3 and £ ; A,- — 1 by definition. These terms 
are used to define the three observables as 



S= -(A2 + A3), 
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Sphericity, Eq. (3), and transverse sphericity, Eq. (4), mea- 
sure the total transverse momentum with respect to the spheric- 
ity axis defined by the four- momenta used for the event shape 
measurement (specifically, the first eigenvector). The allowed 
range of S values is < S < 1, but due to the inclusion of the 
smallest eigenvalue, A3, the typical maximum achieved ex- 
perimentally is S ~ 0.8. Conversely, the transverse spheric- 
ity is constructed using the two largest eigenvalues, and the 
typical range coincides with the allowed range, < S± < 1 . 
Aplanarity (Eq. 5) measures the amount of transverse mo- 
mentum in or out of the plane formed by the two leading jets 
via only the smallest eigenvalue of M xyz , A3, with allowed 
values < A < 1/2. Typical measured values lie between 
< A < 0.3, with values near zero indicating relatively pla- 
nar events. The transverse thrust, T±, and its minor compo- 
nent, T m j_, define a so-called thrust axis for the event, with 



respect to which, the total transverse momentum of the jets 
used in the measurement is minimized. These quantities are 
defined as 
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where 7j_ is translated into Tj_ in order to maintain a com- 
mon event shape definition in which a large value indicates 
a departure from a two-body system. The unit vector fi± de- 
fines the thrust axis of the event. The so-called event plane 
is defined by fi± and the beam direction and allows a mea- 
surement of T m j_. The variable T„ h ± quantifies the sum of 
all transverse momenta p Tl out of the event plane, where the 
sum again runs over each jet i considered in the final state. 
The allowed values for T± span the range < T± < 1 /3 due 
to the range over which both 7j_ and T m i _ may fall, < 7j_, 
T m , X < 2/3. 

Event shapes constructed using hadronic jets in this way 
offer several advantages over explicit cross-section calcula- 
tions for inclusive and multi-jet production. Event shapes 
may be defined as normalized ratios of hadronic final state 
observables, thus reducing the sensitivity to experimental 
uncertainties. Various choices of event shape quantities can 
also lead to enhanced or suppressed sensitivity to differ- 
ent components of the fundamental physical processes in- 
volved [11]. The effect of the underlying event and parton 
shower can be reduced by focusing only on the leading jets. 
The choice of renormalization and factorization scales used 
in calculating the LO and NLO cross-sections may be less 
important when considering ratios of quantities. Systematic 
uncertainties, such as the jet energy scale and detector ef- 
fects, are mitigated by examining the normalized shapes as 
opposed to absolute cross-sections. 

3 The ATLAS Detector 

The ATLAS detector [19, 20] provides nearly full solid an- 
gle coverage around the collision point 1 with an inner track- 
ing system covering |tj | < 2.5, electromagnetic and hadronic 
calorimeters covering |tj| < 4.9, and a muon spectrometer 
covering |tj | < 2.7. Of the multiple ATLAS subsystems, the 
most relevant to this analysis are the inner tracking detector 
(ID) [21], the barrel and end-cap calorimeters [22, 23] and 
the trigger [24] . 



1 ATLAS uses a right-handed coordinate system with its origin at the 
nominal interaction point (IP) in the centre of the detector and the z- 
axis along the beam pipe. The x-axis points from the IP to the centre 
of the LHC ring, and the y axis points upward. Cylindrical coordinates 
(r, <j>) are used in the transverse plane, <j) being the azimuthal angle 
around the beam pipe. The pseudorapidity is defined in terms of the 
polar angle 6 as TJ = — lntan(0 /2). 
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The ID is comprised of a pixel tracker closest to the 
beamline, a microstrip silicon tracker, and lastly a straw -tube 
transition radiation tracker at the largest radii. These systems 
are layered radially upon each other in the central region. A 
thin solenoid surrounding the tracker provides an axial 2T 
field enabling measurement of charged particle momenta. 

The calorimeter is built of multiple sub-detectors with 
several different designs spanning the pseudorapidity range 
up to 1 77 1 < 4.9. The measurements of event shapes are pre- 
dominantly performed using data from the central calorime- 
ters, comprised of the liquid argon LAr barrel electromag- 
netic calorimeter ( 1 77 1 < 1 .475) and the Tile hadronic calorime- 
ter (| 77 1 < 1.7). Three additional calorimeter subsystems are 
located in the forward regions of the detector: the LAr elec- 
tromagnetic end-cap calorimeters, the LAr hadronic end-cap 
calorimeter, and the forward calorimeter comprised of sepa- 
rate electromagnetic and hadronic components. 

The precision and accuracy of energy measurements made 
by the calorimeter system is integral to this analysis and the 
procedures to establish such measurements are described in 
Ref. [25]. The baseline electromagnetic (EM) energy scale 
of the calorimeters derives from the calibration of the elec- 
tronics signal for the energy deposited by electromagnetic 
showers. The hadronic calorimeter has been calibrated with 
electrons and muons in beam tests and the energy scale has 
been validated using muons produced by cosmic rays with 
the detector in situ in the experimental hall [23]. The invari- 
ant mass of the Z boson in Z — > ee events measured in situ in 
the same data-taking period is used to adjust the calibration 
for the EM calorimeters. 

Dedicated trigger and data acquisition systems are re- 
sponsible for the online event selection which is performed 
in three stages: Level 1, Level 2, and the Event Filter. Level 1 
utilizes information from the calorimeter and muon systems 
using hardware-based algorithms. Level 2 and the Event Fil- 
ter are collectively referred to as the High Level Trigger and 
utilize software algorithms running on large farms of com- 
mercial processors. The measurements presented in this pa- 
per rely primarily on the hardware-based Level 1 calorime- 
ter trigger. At this level, coarse calorimeter information is 
used to reconstruct jets in the trigger system with a square 
sliding-window algorithm in 77 — (j) space. 



4 Data Samples and Event Selection 

4.1 Data Sample and Event Selection 

The data used for the analysis of event shapes represent the 
entire 2010 dataset collected at y/s = 7 TeV and correspond 
to an integrated luminosity of / «ifdf = 35. 0± 1.1 pb _1 [26]. 

A sample of events containing high-£>T jets is selected 
via a Level 1 inclusive single jet trigger with a nominal trans- 
verse energy threshold of 95 GeV at the EM scale. The of- 
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Fig. 1 Inclusive jet trigger efficiency as a function of jHj 2 evaluated 
offline. The efficiency to select events with j//t.2 > 250 GeV using 
this trigger is greater than 99.8%. 



fline selection requires two leading jets with a mean trans- 
verse momentum |.Ht,2 > 250 GeV and that the rapidity of 
each leading jet satisfy |y| < 1.0. Subleading jets yield non- 
zero values of each event shape variable and must have pj > 
30 GeV and be within |y| < 1 .5 in order to be used in the cal- 
culations. This choice of event selection is partially driven 
by the trigger threshold which is fully efficient only at an of- 
fline jet transverse momentum of pj > 250 GeV. High mo- 
mentum jets are also less susceptible to the impact of mul- 
tiple proton-proton interactions (pile-up) and have a smaller 
jet energy scale relative uncertainty. The use of \Hj ^ has 
been shown to be significantly more stable against higher- 
order corrections to the jet cross section [27]. The inclusive 
single jet trigger efficiency is evaluated with respect to the 
offline \Hj2 selection, as shown in Fig. 1, and is on aver- 
age greater than 99.8%, thereby removing the need for trig- 
ger efficiency corrections. This determination is made in situ 
using a trigger selection with a threshold of 30 GeV at the 
EM scale. 

The presence of pile-up in these data has the potential 
to impact both event selection and jet reconstruction. Ex- 
perimentally, reconstruction of primary vertices using tracks 
measured in the ID provides a measure of the multiplicity of 
such additional interactions on an event-by-event basis. The 
2010 data contain an average number of primary vertices, 
AW, of approximately (iVpv) = 3, with a tail extending to 
N-p\r > 10. The vertex with the highest total squared track 
momentum, £ p\ track , is assumed to be the vertex at which 
the hard scattering that triggered the event occurred. 

Two primary effects are expected from pile-up: augmen- 
tation of the jet energy scale for jets produced in the hard 
scattering and pile-up jets produced directly by the addi- 
tional pp collisions within the same bunch crossing. The 
consequence of the former is typically an offset to the energy 
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scale which is corrected as described below. The presence 
and impact of pile-up jets on the event shape measurements 
is discussed in more detail in Section 5.2. 

4.2 Jet Reconstruction and Calibration 

Jets reconstructed with the anti-fc, algorithm [28, 29] are 
used for the event shape measurements presented here. This 
algorithm yields regular, approximately circular jets whose 
boundaries are well described by the nominal jet radius. A 
jet radius of R = 0.6 is used here, as in many Standard Model 
jet-physics measurements in ATLAS, compared to the smaller 
radius R = 0.4 jets used for a variety of new physics searches 
and top-quark measurements. This choice is made in order 
to minimize jet-by-jet corrections due to higher-order emis- 
sions and to maximize the reconstruction efficiency since 
the additional small uncertainties that accompany the larger 
radius do not significantly impact the measurement. Anti-fc, 
jets have been shown to be less susceptible than other jet 
algorithms to systematic effects such as pile-up and close- 
by jet activity. The inputs to the jet algorithm are topologi- 
cal energy clusters at the EM scale [30]. Following an aver- 
age offset correction derived in situ to account for noise and 
contributions due to pile-up, an rj- and px-dependent jet en- 
ergy correction referred to as the EM+JES correction [25] 
is applied to all jets to compensate for energy loss in the 
calorimeters, detector geometry and other effects. Jet quality 
criteria such as the timing of calorimeter cell signals, the EM 
energy fraction, and pulse shape information are used to se- 
lect only those events with well measured jets. These criteria 
are designed to remove events that are likely to have con- 
tamination due to beam-related backgrounds, cosmic rays, 
or detector defects. In order to further reduce non-collision 
backgrounds, each event must contain at least one primary 
vertex consisting of at least five tracks with transverse mo- 
menta pj ack > 150 MeV. The MC simulation is reweighted 
in order to match the primary vertex multiplicity observed 
in the data. 

4.3 Monte Carlo Simulation 

Dijet and multi-jet events are generated using two approaches. 
The first uses direct perturbative calculation of the tree-level 
matrix elements in powers of the strong coupling constant, 
as- The matrix elements are evaluated at LO in as for each 
relevant partonic subprocess. This is a so-called "multi-leg" 
method. The second approach implements a sampling of the 
phase space available for gluon emission with some suitable 
approximations. The latter uses LO perturbative calculations 
of matrix elements for 2 — > 2 processes and relies on the 
parton shower implementation to produce the equivalent of 
multi-parton final states. 



The multi-leg technique is used by Alpgen [31]. In the 
analysis presented here, Alpgen 2.13 is used with up to 
six final-state partons. ALPGEN is interfaced to both Her- 
WIG 6.510 [32] to provide the parton shower and hadroniza- 
tion model, and to Jimmy 4.31 [33] for the underlying event 
model. The CTEQ6L1 LO [34] parton distribution function 
(PDF) with LO as is used for ALPGEN. 

The parton shower simulation programs Pythia [35] 
and HERWIG++ [36] both implement the second approach 
for QCD jet production and rely on the parton shower to gen- 
erate multi-jet final states. Pythia 6.423 with the Perugia 
2010 tune [37] and HERWIG++ 2.4.2 are used to compare to 
the data; these provide shower models that are pj ordered 
and angular ordered, respectively. LO PDFs are taken from 
the MRST2007 LO* [38, 39] PDF for HERWIG++, and from 
the CTEQ6L1 LO [34] PDF in PYTHIA. 

The MC programs used for comparison to the measure- 
ments of event shapes are chosen in part for their ability to 
describe other ATLAS jet-based measurements. The multi- 
jet cross section measurements [40], which constitute the 
same final states as those probed in this analysis, exhibit 
very good agreement with the predictions from ALPGEN. 
HERWIG++ not only exhibits good agreement with individ- 
ual jet shape measurements [41] but is also tuned to yield 
good agreement with event shape measurements from the 
Large Electron Positron (LEP) collider experiments. Finally, 
the Perugia 2010 tune of PYTHIA also shows good agree- 
ment with the ATLAS jet shape measurements and has been 
tuned using the theoretical input from higher-order calcula- 
tions of event shapes presented in Ref. [11]. These three MC 
programs thus provide well motivated predictions for the fi- 
nal state observables measured via event shapes. 

Events generated by these MC programs are passed through 
a full simulation [42] of the ATLAS detector and trigger 
based on GEANT4 [43] and processed in the same way as 
collision data. The Quark-Gluon String Precompound (QGSP) 
model [44] is used for high energy inelastic scattering of 
hadrons by nuclei, and the Bertini cascade model [45] is 
used to describe the interactions of hadrons with the nu- 
clear medium. Alternative GEANT4 physics lists that spec- 
ify particle and process definitions, using a combination of 
the FRITIOF [46] and Bertini models and QGSP without 
Bertini, are used as part of the studies to understand the un- 
certainties on the jet energy scale. 

5 Results and Systematic Uncertainties 

The event shape measurements using jets presented here are 
corrected to particle-level, after accounting for detector effi- 
ciencies and instrumental effects. Particle-level jets are con- 
structed from all final state particles from the MC simula- 
tion with lifetimes longer than 10 ps. Direct comparisons 
can thus be made between the results presented here and 
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MC generator data after parton shower and hadronization. 
The dependence of these observables on iffr,2 is also eval- 
uated. This allows both the investigation of the results and a 
clearer picture of regions of phase space in the measurement 
where MC predictions may not fully describe the data. 

5.1 Accounting for Detector Effects 

In order to compare the predictions of MC event generators 
with the measurements, several effects must be accounted 
for. Efficiency loss due to detector coverage and resolution, 
small detector biases, and other effects, may all affect the 
measured value of an event shape variable. In order to ac- 
count for such effects, the MC and detector simulation are 
used to estimate their impact. MC events after full detec- 
tor simulation are used to derive bin-by-bin corrections that 
are applied to the detector-level measurements of each event 
shape variable to obtain the unfolded, particle-level result to 
which the MC simulated events after parton shower may be 
compared. These corrections are typically smaller than 10%. 
The bin sizes were chosen to be approximately commensu- 
rate with resolution, with individual bin purities significantly 
greater than 50%. 

The primary MC generator used for evaluating the cor- 
rections is Alpgen, since the detector-level distributions 
are well described and Alpgen models the multi-jet cross- 
section well [40]. As a cross-check of the method, the cor- 
rections evaluated with ALPGEN are compared to those ob- 
tained from HERWIG++ and Pythia. For \ny23, A, and S, 
this component of the uncertainty is approximately 2%-8%, 
which is smaller than both that due to the overall JES sys- 
tematic uncertainty discussed below and the finite sample 
size with which the correction factors are determined. How- 
ever, for the thrust event shape variables, in particular for 
T m x, and S±, the generator dependence of these corrections 
is approximately 10% for the majority of the range of those 
measurements. 

5.2 Systematic Uncertainties 

Multiple effects are present in the measurement of event 
shapes due to the inclusive nature of these observables. These 
include the uncertainty due to the jet energy scale, the ef- 
fects of multiple pp interactions, the finite resolution, and 
the fiducial range of the detector. All of these effects are 
evaluated and accounted for in the measurement. The domi- 
nant uncertainties are the jet energy scale and generator de- 
pendence of the corrections in regions of high statistical pre- 
cision. 

The uncertainty on the JES established by the jet calibra- 
tion procedure [25] influences the final event shape measure- 
ment via both the thresholds used to select events and the 



momenta used to calculate the event shape observables. This 
uncertainty is primarily established by the measurement of 
the single hadron response using test beam data, but is also 
verified in situ during 2010. For jets used in these measure- 
ments the typical JES uncertainties are 2.3%-3.0%. The im- 
pact of this source of systematic uncertainty is reduced by 
the explicit use of ratios of jet momenta for several observ- 
ables, although the jet yield can still vary for a given event 
due to these selections. Variations of the individual jet mo- 
menta are performed within the systematic uncertainties of 
the JES measurement. For nearly all measured event shape 
variables, the overall JES uncertainty has the largest impact 
apart from statistical precision. Most observables have an 
approximate 5% uncertainty due to the JES, with A and T± 
being impacted by up to 10% in the steeply falling tails of 
the distributions. 

Additional jets present in the event due to pile-up may 
also alter the observed event shape. This may be of partic- 
ular importance for those measurements that are explicitly 
dependent upon the jet multiplicity, such as those computed 
from the event transverse momentum tensor. The impact of 
pile-up on the jet energy is accounted for by the energy scale 
corrections and uncertainty discussed above, where the pile- 
up contribution becomes negligible for pi > 250 GeV. 

A crucial tool in the identification of jets from pile-up 
is the jet-vertex fraction, or JVF [47]. This discriminant es- 
timates the contribution of pile-up to a single jet by mea- 
suring the fraction of charged particle momentum in the jet 
that originates in the hard scatter. The rate of additional jets 
from pile-up in events with between five and eight primary 
vertices exhibits an increase of a factor of two compared to 
events with two primary vertices. Because the overall frac- 
tion of events with greater than five reconstructed primary 
vertices is approximately two percent, and these jets tend 
to have a much softer pj spectrum, the impact due to addi- 
tional jets is significantly smaller than other systematic and 
statistical uncertainties for all measurements. 

To further establish the systematic uncertainty incurred 
by pile-up, comparisons are made between the observed detector- 
level distributions in events with and without additional re- 
constructed vertices. A slight variation of a few percent at 
low ln>'23 is observed as well as a relative 10% increase in 
the fraction of events at higher T±. Similar observations are 
made by evaluating the impact of the JES uncertainty. Fur- 
thermore, each event shape is measured as a function of the 
JVF of the third jet in the event to directly test the effect of 
pile-up on the final state observables. Variations of less than 
3% are observed when requiring that jets contain a high frac- 
tion of associated track momentum originating in the identi- 
fied hard-scatter vertex. This effect is taken into account in 
the systematic uncertainty in the final result. 



6 



5.3 Event Shape Distributions 

The normalized distributions of the third-jet resolution pa- 
rameter and aplanarity are shown in Figs. 2(a) and (b). In the 
case of y23, where the primary sensitivity is to the descrip- 
tion of the momentum of the third jet, Pythia provides the 
most accurate description of the data, whereas HERWIG++ 
exhibits slightly better agreement than Alpgen. Although 
Alpgen provides exact tree-level matrix element calcula- 
tions for up to six jets, it overestimates the fractions of events 
in the range lny23 < —5. This is qualitatively expected be- 
cause Alpgen's more precise calculation of the high jet 
multiplicity states is primarily concerned with jets near the 
hard scale of the event. In this region, the impact of the 
leading-logarithm resummation in Pythia and HERWIG++ 
plays a significant role. As a result, PYTHIA and HERWIG++ 
both describe the data more accurately than ALPGEN for this 
event shape variable, particularly at small values of In V23. 

Aplanarity measures the sum of the transverse momenta 
out of the event plane defined primarily by the two hardest 
jets. The deviation of the MC prediction from the data is 
significant for HERWIG++, with some small differences ob- 
served with respect to ALPGEN as well. The measurements 
consistently support more highly aplanar events than pre- 
dicted by HERWIG++, with the majority of the distribution 
observed to be significantly different from the MC predic- 
tion. The agreement with PYTHIA is good across the full 
distribution. These results suggest that the event shape is 
more accurately described by the exact multi-jet prediction 
provided by the multi-leg matrix element generator (such as 
Alpgen) and the model provided by Pythia. 

The statement that the topological distribution is bet- 
ter described by Alpgen and Pythia than by Herwig++ 
is supported by the measurement of the transverse thrust, 
Tj_. Fig. 2(c) exhibits the same behavior as observed in the 
aplanarity: HERWIG++ predicts fewer than observed highly 
isotropic events at large Tj_. Throughout the distribution, ALP- 
GEN and Pythia both predict the measured thrust well. The 
minor component of the thrust, T m ±, or the out-of -plane 
thrust magnitude, shown in Fig. 2(d), does not exhibit as 
large a difference as observed in the aplanarity. A slight 
overestimation by PYTHIA is observed for intermediate val- 
ues of the thrust minor component, 0.25 < T m j_ < 0.40. 

Lastly, the sphericity and transverse sphericity distribu- 
tions shown in Figs. 2(e) and (f) exhibit differences between 
all three generators and the data. The construction of the 
transverse sphericity as a ratio of eigenvalues of the mo- 
mentum tensor of the event leads to a slightly improved de- 
scription. In both cases, PYTHIA provides the best descrip- 
tion of the data, and the data are better described by ALP- 
GEN than HERWIG++. In particular, HERWIG++ underesti- 
mates the number of highly spherical events in the range 
0.40 < S < 0.72. 



5.4 Dependence on i^fr,2 

The measurement of the distribution of event shape observ- 
ables allows for a detailed comparison with MC predictions 
for a large range of kinematic phase space defined by j//t,2 
in multi-jet events. It is informative to evaluate the explicit 
dependence of these shapes on the kinematic properties of 
the event in order to determine potential differences in the 
modeling of the data for different jet momentum ranges. The 
evolution of each event shape variable with 5//T.2 exhibits 
a similar trend as that expected by the running of as which 
leads to a reduction in extra gluon radiation and thus a re- 
duction in the value of each event shape. 

Figure 3 depicts the dependence of the mean of each 
event shape variable on |Ht,2 as it approaches 1 TeV. In 
all cases, a general trend is observed in which the mean de- 
creases as i/fr,2 increases. This can be understood in terms 
of how the dynamics of the 2 — > 2 process evolves with en- 
ergy. As the energy in the leading jets increases, the di-jet 
structure dominates because of kinematics and because as 
decreases as j^t.2 increases, causing the dominant NLO 
corrections which generate higher relative momentum gluon 
emission to decrease as well. 

The variation as a function of ji£r,2 is the largest for 
(In 723), which indicates a change of j23 of nearly a factor of 
five between j#t,2 = 300 GeV and 800 GeV. Nonetheless, 
the agreement between the MC prediction and the observed 
event shape variable dependence is good for all generators. 
The observation made above that too few highly aplanar 
events are present in the predictions from HERWIG++ is 
again observed in the evolution of (A) with the momentum 
scale of the event. The agreement among the mean values 
measured and the MC predictions improves for j//t,2 > 
500 GeV, although the systematic uncertainties are larger 
and the statistical power of the measurement is reduced. 
Similarly, the evolution of (t_l) with j//t,2 in Fig. 3(c) is 
underestimated by HERWIG++. ALPGEN and PYTHIA con- 
sistently predict the mean value of T± and its evolution with 
jH-y.2 more accurately, whereas all of the generators de- 
scribe (T m jl) vs. \Hi 2 well. Finally, the sphericity and the 
transverse sphericity (Figs. 3(e) and (f)) are both measured 
to be approximately 10% larger than predicted by the three 
MC programs at low j.Ht,2, whereas the agreement again 
improves for higher values. This difference is driven by the 
underestimate of highly spherical events observed in Figs. 2(e) 
and (f) which decreases as the average sphericity decreases 
as a function of \Hj2- 

6 Summary 

Six event shape observables are measured with jets in proton- 
proton collisions at y/s = 7 TeV with the ATLAS detector 
with a data sample of 35 pb . 
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Fig. 3 Mean value of each event shape variable as a function of \Hj,2- Comparisons are made between the MC generators HERWIG++, ALPGEN 
and Pythia. 
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Measurements are performed up to an event Ht.2 of 2 TeV 
and are compared with different Monte Carlo event gener- 
ators. Overall shape comparisons are made with these MC 
programs, as well as the kinematic evolution of the mean 
value of each event shape variable with \Hj2- Reasonable 
agreement is observed in most kinematic and topological re- 
gions. The measurements suggest that the modeling of the 
data by Pythia (Perugia 2010) and Alpgen are more ac- 
curate than that by HERWIG++, in particular for the apla- 
narity, A, and transverse thrust, Tj_. The good description 
provided by Alpgen (+Herwig/Jimmy) of the multi-jet 
cross-section [40] is reflected as well in these measurements, 
although the description provided by Pythia tends to model 
the data more accurately. PYTHIA predicts a slightly higher 
mean T m j_ at low ji/j,2 than observed in the data whereas 
HERWIG++ predicts a slightly lower mean thrust. The sys- 
tematic uncertainties in the measurement of S and S± are 
found to be relatively small. The observation that the mea- 
sured mean value of each event shape variable decreases 
with jHj 2 is consistent with the trend expected from the 
running of as and is generally well modeled by the MC sim- 
ulations. 

Comparisons of these results to previously published LHC 
measurements of hadronic event shapes [7] indicate that the 
slight overestimate by PYTHIA of events with thrust minor 
component, T m ±, in the range 0.25 < T m j_ < 0.40 is ob- 
served in each case. However, the good agreement observed 
in this study between data and both PYTHIA and ALPGEN 
for the thrust, T± is not seen in Ref. [7]. The different tunes 
of HERWIG++ and PYTHIA, as well as the different under- 
lying event and hadronization models interfaced to ALP- 
GEN in the two measurements, may account for these differ- 
ences. Furthermore, the systematic uncertainties associated 
with these measurements are in many cases similar to the 
differences between the generators. 

ATLAS measurements of the jet shape [41], jet fragmen- 
tation function [48], and multi-jet cross-section also pro- 
vide additional insight into the results shown here. PYTHIA 
and HERWIG++ 2.4.2 both provide a reasonable descrip- 
tion of the fragmentation function of high momentum jets 
in the data, whereas only the former models the jet shape 
accurately. ALPGEN, on the other hand, yields a fairly ac- 
curate description of the three-jet to two-jet cross-section 
ratio [40], although it produces jet shapes that are signifi- 
cantly narrower than those measured here. 

These results, supported by other ATLAS measurements 
of the hadronic final state, reinforce the importance of the 
leading-logarithm resummation in parton shower MC event 
generators. They also demonstrate the ability of leading or- 
der MC to provide a reasonable description of multi-jet event 
shapes. 
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